Synthetic Observations of Carbon Lines of Turbulent Flows 
in Diffuse Multiphase Interstellar Medium 
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ABSTRACT 

We examine observational characteristics of multi-phase turbulent flows in the diffuse inter- 
stellar medium (ISM) using a synthetic radiation field of atomic and molecular lines. We consider 
the multi-phase ISM which is formed by thermal instability under the irradiation of UV photons 
with moderate visual extinction Ay ~ 1. Radiation field maps of C"*", C", and CO line emissions 
were generated by calculating the non-local thermodynamic equilibrium (nonLTE) level popula- 
tions from the results of high resolution hydrodynamic simulations of diffuse ISM models. By 
analyzing synthetic radiation field of carbon lines of [C II] 158 /J,m, [C I] ^P2 Pi (809 GHz), 
^Pi Pq (492 GHz), and CO rotational transitions, we found a high ratio between the lines of 
high- and low-excitation energies in the diffuse multi-phase interstellar medium. This shows that 
simultaneous observations of the lines of warm- and cold-gas tracers will be useful in examining 
the thermal structure, and hence the origin of diffuse interstellar clouds. 

Subject headings: ISM: clouds - ISM: molecules - radio lines: ISM 



1. Introduction 

Recent studies of the multi-phase interstellar 
medium (ISM) have demonstrated that the ori- 
gin of small scale interstellar turbulence might 
be attributed to the thermally unstable na t ure of 
the ISM itself (iKoyama fc Inutsukal 12001 l2006t 



Kritsuk fc Norma n' '2002a'! bl: lAudit fc Henneb'eiij 

2005. : ,Heitsch et al.„2005;,Inutsuka. Kovama fc Inoud 



20051 : IVazauez-Semadeni et all |2006[) . Until re- 
cently, most theoretical work on interstellar tur- 



bulence has assumed that the ISM is isothermal 
(see for recent reviews, Mac Low fc Klessen 11200-5 



Elmegreen fc Scald |2004 , and references therein). 
This approximation has been imposed mainly be- 
cause quiescent gas evolves almost isothermally 
at lOK as a result of the thermal balance be- 
tween radiative heating and dust emission cool- 
ing where external UV radiation i s well-shielded 
( Hattori. Nakano fc Havashi 1969f ). and partly 
because of limited computational ability. De- 
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tailed analyses of the one-dimensional collapse 
due to thermal instabilit y can be found in, for 
example, .Hcnncbcllc fcj ^ault I ([l999, 2000) and 
Kovama fc Inutsukal (feoool) . The situation has 



changed since Kovama &: Inutsukal (|2002l ) first 
pointed out that the thermally unstable nature 
of multi-phase ISM is able to maintain the tur- 
bulent flows. Since then, various attempts have 
been made to examine the properties of the multi- 
phase turbulent ISM dKqyama fc Inutsuka 20021: 



Kritsuk &: NormanI l2002allbl : [Audit fc Hennebelld 
2005) . Kovama fc Inutsukal (|20 02) performed two- 
dimensional simulation of a supernova-induced 
shock propagating into warm (~ 6000K) ISM. 
They demonstrated that broad emission linewidths 
can be generated in the post-shock layer, where 
chaotic structures are formed by thermal insta- 
bility. In their calculation, tiny, cold and dense 
clumps were seen to move in the surrounding 
warm diffuse medium subsonically with respect 
to the warm medium, but supersonically to the 
sound speed in the cold clumps. In this article, 
we refer to this model as the "two-phase model" 
where two thermal states, i.e., cold clumps and 
warm medium, are generated by thermal instabil- 
ity. Similar conclusions have a lso been reached by 
other a uthors. Among them, iKritsuk fc Norman 
(|2002al lbh ])erformed three dimensional simula- 
tions of a multi-phase medium and examined the 
evolution of kinetic energy of t he re sulting tur- 



bulence. lAudit fc Hennebelld (j2005[ ) simulated 
a two-dimensional converging flow, where the 
post-shock thermally unstable gas is continuously 
provided by the two shock waves that confine 
the post-shock gas. They s howed that complex 
struct ures similar to those of iKovama fc Inutsuka 
(|2002l ) are formed, and discussed their evolution 
by comparing cooling length, fragment scale of the 
cold neutral medium (CNM) , Field l ength, and the 
typica l size of the shocked layer. Hcitsch et al 



(|2005l ) and IVazquez-Semadeni et al.i ((2006) have 
followed this line of work and discussed the ob- 
servational signature as well. While the calcula- 
tions above focus ma. inly on neutral atomic gas, 
Pavlovski etldl (|2005h examined molecular gas of 
much higher density (nn — 10^ cm"'^). They con- 
cluded that the gas is kept isothermal due to the 
fast reformation of the coolant molecules behind 
shock waves in the dense molecular clouds. 

In this article, we concentrate on more dif- 



fuse gas (w 1 cm~'^), which is warmed by mild 
UV photoelectric heating, and focus on observa- 
tional characteristics of the two-phase medium. 
Emission-line diagnosis of atomic and molecular 
radiation in millimeter and submillimeter bands 
provides a powerful tool for o bservat ional studj^of 
interstellar turbul ence (e.g., IPadoan et al.l I2OOII : 
Bergin et al. 20041 ). In addition to existing tele- 
scopes equipped with submillimeter receivers, such 
as SubMillimeter Array (SMA) and Atacama Sub- 
millimeter TElescope (ASTE), the highly sen- 
sitive instruments of Atacama Large Millime- 
ter/submillimeter Array (ALMA) will be available 
before long. We thus examine how a characteristic 
feature of the two-phase turbulence model should 
appear in observational quantities. We try to 
find and propose unambiguous and robust obser- 
vational indicators for the two-phase turbulence 
model. For this purpose, we perform hydrody- 
namical simulations of two-phase turbulent flows, 
and then calculate line emission of atomic carbons 
and CO molecules. Based on the synthetic emis- 
sion maps, our results are more suitable for direct 
comparison with observation than the density or 
temperature distributions, which are provided by 
previous hydrodynamical simulations. 

The organization of this article is as follows: In 
§2 we describe our model and basic equations. In 
§3 we display our synthetic radiation maps of the 
two-phase model as well as the conventional one- 
phase (isothermal) turbulent flow model. We also 
show the results of the analyses of our simulated 
emission maps of C"*", C", and CO lines. In §4 
we discuss the limits of and uncertainty about our 
calculation and discuss its relevance to current ob- 
servational results. We summarize this article in 
the final section (§5). 

2. Models and Calculations 

We consider diffuse clouds irradiated by mild 
UV background radiation of Gq = 1 (where Go 
is the typical interstellar UV strength in Habing 
unit), which is attenuated with average visual ex- 
tinction Ay < 0.1. We perform hydrodynamical 
simulations first, and then calculate line-cmissivity 
distributions using the simulation data. In the fol- 
lowing, we present equations and assumptions for 
hydrodynamical and radiation calculations in this 
order. 
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2.1. Hydrodynamic Equations 

We solve an ordinary set of hydrodynamical 
equations: 



dp 

dt 
dv 

m 



+ V-{pv) = 0, 
+ {v ■ V)u = 



^ + V • {pev) +PV-V - - ( ^ I A(T) 



(la) 
(lb) 

V • {KVT) (Ic) 



VP 
P 

P 



where p is the density, P is the pressure, e = 
P/(7 — 1) is the specific internal energy with 
the adiabatic index 7 = 5/3, T and A repre- 
sent the heating and cooling rates, respectively, 
and V is the velocity of a fluid element. The last 
term of the third equation (jlcp represents thermal 
conduction. The thermal conduction coefficient 
K = 2.5 X lO^r^/^ ergs cm^i K^^ sec'^ is taken 
from iParken (|l953r ). which is valid for tempera- 
ture T < 4.47 X lO^'K. For cooling and heating 
rates, we empl oy the fitting formulae sho wn by 
equation (8) of Kovama fc Inutsukal ( 20061 ). The 
formulae are de rived by approximately de scrib- 
ing the result of Kovama fc Inutsuka ( 2000l ). who 
performed non-equilibrium thermal calculations of 
one-dimensional shock waves. Their formulae con- 
tain dust-grain photoelectric heating, cosmic ray 
and X-ray ionization heating, H2 formation and 
destruction, and line cooling by atomic and molec- 
ular emission. In equation (jlbp we do not intro- 
duce the physical viscosity, which affects the sat- 
uration level of turbulent motions, because we are 
more interested in characteristic features in radi- 
ation field rather than detailed arguments on the 
mechanism for maintaining the turbulent motion. 
Detailed systematic analyses of the saturation am- 
pli tude of two-phase tu r bulen t motions are found 



Kovama fc Inutsukal (|2006[ ). We can safely ig- 



nore back reactions by cooling radiation to hydro- 
dynamical evolutions either via radiation pressure 
or radiation force. 

Our two-dimensional simulations employ a 
scheme based on the second-order Godunov 
method. In order to correctly resolve thermal 
conduction-induced structures, the Field length 
Ap = (i^T/n^A)^/^, which is the characteristic 
length scale of thermal conduction, must be re- 
solved by more than three cells (the so-called Field 



condition, after ama fc Inutsukal |2004| ). Our 
calculations satisfy this condition with the num- 
ber of Eulerian grids 2048^ in the simulation box 
of 2.4x2.4 pc^. 

We introduce a supersonic velocity field to the 
medium of initially uniform density and temper- 
ature, and calculate its evolution without an ex- 
ternal driving force (i.e., so-called decaying turbu- 
lence). The initial supersonic velocity field has an 
average Mach number of (Ai) « 10, and power 
spectrum P{k) (x in the range of 1 < fc/27r = 
1/A < 4, in units of the reciprocal of the size of 
the simulation box A. The purpose of this paper is 
to simulate the radio band observation of certain 
samples of turbulent clouds. Thus we set up the 
initial conditions rather arbitrarily and calculate 
the dynamical evolution of the relaxation process 
from these highly turbulent flows. For the actual 
simulation of the observation described below, we 
choose snapshots of the time evolution that would 
mimic the observed velocity proflles (i.e., position- 
velocity diagram) of diffuse clouds. 

2.2. Method of Radiation Calculation 

We solve non-local thermodynamic equilibrium 
(nonLTE) level populations for relevant lines of 
atomic and molecular carbon, 0+, C*^, and CO. 
Each level population is determined by solving the 
equations of detailed balance for transitions as a 
function of local density and temperature taken 
from the hydrodynamical calculations. In our 
model of mildly UV-irradiated diffuse medium, 
in addition to the Lya line of atomic hydrogen 
in the UV band, various millimeter and submil- 
limeter lines contribute to the cooling. Among 
the line-cooling processes, the contribution of 
[C II] fine structure cooling is dominant in a 
wide densit y range of diffuse regime s (see Fig- 
ure la of Kovama fc Inutsuka 2000l ). Toward 
higher density, the role of C° and CO coohng 
gradually increases. Around the column density 



10 cm where background UV radia- 



tion of Go = 1 begins to be shielded, the domi- 
nant coolant changes from ionic to molecular car- 
bon. Our ISM model, which includes heating and 
cooling under Nh — 10^"cm"^, corresponds to 
the atomic/molecular transition region, and thus 
ionic, neutral and molecular carbon species should 
exist simultaneously. Therefore we examine three 
species of carbon, i.e., C+, C°, and CO, in order to 
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probe various thermal states in multi-phase ISM. 
We calculate C II and C I fine structure transi- 
tions, whose excitation energies are 92K([C II]), 
24K([C If Pi -3 Po), 39K([C lfP2 Pi), and 
63K([C I]3P2 -3 Po), as well as CO(J-(J - 1)) 
rotational transitions, with 5.5 J K, where J is the 
rotational quantum number. 

For [C II] 158Mm lin e, wc use a modified cooling 
function bv lKovama fc Inutsuka (|2002. ]) : 



A 



CII 



2.8 X IQ-^^Vt ■ cxp (^-jr^ 
xy(C+)[y(H0)+y(H2)]nH' 



erg sec ^cm 



(2) 



where y{X) = n{X)/nii represents the fractional 
abundance of the line-emitting particle X. Both 
atomic and molecular hydrogen work as collision 
partners for exciting C"*" . In the regime where C II 
cooling dominates, we can safely replace the term 
for collision partners with y(H°) -I- y(il2) = y{ii'^)- 
The emissivities per unit volume by [C I] and 
CO lines are calculated by solving the nonLTE 
level population. The emissivity of species X by 
transition from level i to j per unit volume (cool- 
ing rate) is given by 



^x,2j = hvijAijn{X,i), 



(3) 



where hvij is the energy difference between lev- 
els i and j, Aij is the Einstein A-coefficient for 
spontaneous emission from energy level i to j, and 
n{X,i) is the population density of species X and 
energy level i. The level populations are obtained 
by solving the equations of statistical balance: 

j' j' 

(4) 

In the above equation, Cij is the collisional tran- 
sition rate per emitting particle, which is written 
with the collisional rate coefficient 7^- as 



(5) 



where n{Y) indicates the density of collision 
partner Y. For neutral carbon atoms, we con- 
sider the triplet in the electronic ground state 
('^Po,3 Pi,'^ P2). Relevant transition coefficients 



collision partners for excitation of atomic carbon 
include atomic hydrogen and electron, and ne- 
glect H2 for only a small contribution of II2 is 
expected where C I abundance and line emission 
have their peaks (the abundance of C° reaches its 
maximum between the surface of UV irradiated 
region and the well-shielded inner region where 
most of carbon is locked in CO molecules). For 
CO molecules, we make use of a simp lified formula 
for th e collisional rate coefficients of iMcKee et al. 
(jl982l ). CO is assumed to be excited by collision 



with atomic and molecular hydrogen. We use the 
same collisional coefficient as atomic hydrogen for 
H2. 

In order to calculate line intensities we need 
to know the ionization degree and the chemi- 
cal abundances for each species. This requires 
non-equilibrium calculation of chemical evolution 
along with hydrodynamical simulation, which is 
computationally very expensive. To avoid this 
difficulty, here we simply assume uniform abun- 
dance distribution, and a prescribed variation in 
terms of temperature (see §4.1|) . For the ion- 
ization degree, we assume an electron fraction 
of y{e) = lO^'^ referr ing to the th e rmal equi- 



librium abundances of IWolfire et al.l ( 19951 ) . In 



the neutral medium we are considering, change 
in the assumed value of y(e) does not signifi- 
cantly alter the results since the collisional tran- 
sition by electrons is not dominant. For carbon 
species, we take y(C+) = lO^'^ and y(C°) — 
y{GO) = 3 X lO"*". These abundances are taken 
from typic al values in the one - dimen sional cal- 
culation of iKovama Sz Inutsukal ( 2000h except for 
C°. The abundance of C'' was not calculated in 



Kovama fc Inutsuk a (2000), and we arbitrarily as- 
sumed the same value as CO for simplicity. The 
fractions of C° and CO are affected by such fac- 
tors as the degree of shielding of ultraviolet radi- 
ation and the abundances of other species. We 
try to derive quantities that are independent of 
the abundance distribution, and defer discussion 
on the effects of chemical evolution in our model 
to later in i |4T] 

Line emissivity of each fluid element has the 
thermally broadened profile 



0(^) 



1 



exp 



(6) 



are taken from iHollenbach fc McKed (|l989f) . The 



where vq denotes the frequency of the line center. 
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and Ai^D = {vq/c) ^^ksT jm denotes the ther- 
mal width at temperature T. Line profiles are ob- 
tained by convolving the thermal broadening with 
the line centroid i^o Doppler-shifted by the bulk 
fluid motion. The bulk fluid velocity is taken from 
results of the hydrodynamical calculations. 

3. Results 

3.1. Line Emission Maps 




\Ji -0.5 0.0 fiS 1.0 



Fig. 1. — (a) Temperature and (b) density dis- 
tributions of the two-phase flow model at time 
t = 2.479 Myr. Numerous tiny dense and cold 
clumps are formed within the simulation box. 

Figure [T] shows a snapshot of temperature and 
density distributions. In the density distribution 
map (Figure [T]d) , numerous tiny dense clumps are 
observed, having been formed as a consequence of 
phase transition in the thermally unstable medium 
after the initial shock heating. Comparing Figures 
[T^ and b, we can easily recognize that the den- 
sity map is just a color-inverted version of the 
temperature map. This is due to the fact that 
this turbulent flow is essentially driven by the 
thermal instability characterized by its approx- 
imate isobaricity. Under almost uniform pres- 
sure, thermally unstable gas forms numerous tiny 
cold and dense clum ps within the surrounding 
warm and diffuse g as ("Kovama & Inutsuka"200^ 
Kritsuk &: NormanI |2002a; Audit & Hcnnebelle 



Figures[2^ and[2lD represent the synthesized 
emissivity maps of [C II] 158 fj,m and C0(1 — 0) 
lines. Since the density is below the critical den- 
sity for LTE ncrit(~ 3000cm"3) throughout the 
whole region, and both lines are assumed to be 
optically thin, line emissivity per unit volume is 
proportional to under the assumption of uni- 
form abundance distribution. Thus the emissivity 
map closely resembles the density distribution one. 
Tiny cold clumps strongly emit radiation in both 
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Fig. 2.— Synthesized (a)[C II] and (b) C0(1 - 0) 
line-emissivity distributions in the two-phase flow 
model a.tt = 2.479 Myr. Intensities of the emission 
lines, which are in proportion to the emissivity in 
the optically thin case, reflect both the density and 
temperature distributions. Radiation is emitted 
mainly from tiny, cold, dense clumps. 



[C II] and C0(1 — 0) lines because of high density. 

In order to extract outstanding characteristics 
of two-phase medium, we also perform the sim- 
ulation of "one-phase" isothermal (lOK) turbu- 
lent flow with the same resolution (2048^ grids). 
This choice of temperature is based on the conven- 
tional assumption in many isothermal turbulence 
models of dense molec ular clouds strongly shielded 
against UV field fe.g.. |Padoan et al.ll200l[) . Figure 
[3] shows the density distributions from a snapshot 
of the one-phase model. Filamentary structures, 
which are commonly observed in simulations of 
isothermal turbulent fiow, appear in Figure [3] as 
well. Sharp edges of filamentary structures im- 
ply that they are formed by shock compression. 
The synthetic emission maps of [C II] I58/im and 
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Fig. 3. — Density distribution in the one-phase 
flow model at t = 5.248 Myr. Woven structures 
with sharp edges are observed, formed through 
strong shock compression by initial supersonic tur- 
bulent motions. 
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C0(1 — 0) lines, which are calculated in a way 
similar to those of the two-phase model, are dis- 
played in Figure H) As in the two-phase model, 
[C II] and C0(1 — 0) emission distributions closely 
trace the density one. For the one-phase model, 
our emissivity maps agree well with previous work 
on, for example, synthetic maps of optically thin 



^^CO line of the isothermal model of lPadoan et al 
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Fig. 4. — Emissivity distributions of (a) [C II] and 
(b) CO(I — 0) lines in the one-phase flow model 
with T = lOK at t = 5.248 Myr. Highly woven 
structures are also seen in line-intensity distribu- 
tions. 



3.2. Analyses of Synthetic Line Intensities 

To extract characteristic properties of the two- 
phase models from our calculations, further anal- 
ysis is required. We simulate realistic observations 
of a three-dimensional cube by taking one of the 
axes of our two-dimensional hydrodynamical cal- 
culations as a line of sight. Defining the x axis as 
the line of sight and the x component of the ve- 
locity field as the velocity along the line of sight, 
we can calculate the intensity of the line emission 
I{L) as 



47r dvx 



(7) 



where A^ is the emission rate per unit volume on 
the L-th line of sight, Ax is the grid size, and Avx 
is the bin size of the velocity channel, respectively. 
Straight summation of A^ is applicable because 
we assume all the lines are optical thin. Hereafter 
we analyze integrated intensity data obtained by 
equation Q. Note A^ includes thermal broaden- 
ing as well as bulk fluid velocity (see eq. H]), and 
therefore I{L) correctly reflects the velocity struc- 
ture inherent in the model as well. 



3.3. Ratio of High- and Lov^r- Temperature 
Tracer Lines 

A distinguishing feature of the two-phase tur- 
bulence model induced by thermal instability is 
the coexistence of the gases of diversely different 
thermal conditions. In this section, we investigate 
characteristics of the two-phase model by simul- 
taneous analysis of two kinds of lines that trace 
either high- or low- temperature medium. In the 
following sections, we examine intensity ratios be- 
tween the [C II] line and C0(1 - 0) Hue, CO rota- 
tional lines, and C I fine-structure lines. 

3.3.1. [C IIJ-CO(J^l-O) line ratio 

Figure [5] displays the h^cii] ~ Aco(i-o) ratio (in 
Figure [5] the ratio of \cii\l ^co{i-o) is multiplied 
by a constant factor 3 x 10"'' just for clarity of the 
figure; therefore, in the following arguments, the 
absolute values of longitudinal axis are in an arbi- 
trary unit, and only the relative positions of the 
points in the diagram have meaning). The points 
taken from the one-phase model fall on a straight 
line of a constant value of A[cii]/Aco(i-o) ^ 10~^, 
while those from the two-phase model make a 
slightly bending curve with much larger values of 
^[Cii]/^co(i-o) ^ 10^ this diagram. This dif- 
ference is interpreted as follows: As we assume 
that the medium is optically thin, and all the den- 
sities in our calculations are below the critical den- 
sity for both [C II] and C0(1 - 0) fines (« 3000 
cm~^), the line ratio is simply described in terms 
of A, 



'[cii] 



A 



[CII] 



y(C+)(az;)[cii](T) 



^co(i-o) Aco(i-o) 2/(C0)(cru)co(j=i)(r)' 

(8) 

In equation ([8]) (cru) [cii] and ((Tw)co(j=i) are de- 
termined by the cross sections for collisional exci- 
tations and are dependent only on temperature. In 
our calculations the abundances of CO and C'^ are 
assumed to be uniform, and then equation ([5]) im- 
plies that the ratio /[cii]/-^co(i-o) is a function of 
temperature. In our one-phase turbulence model 
of which temperature is adopted to be constant 
(=10K), the temperature-dependence of equation 
([5]) drops out and the [C II] intensity is linearly 
proportional to the C0(1 — 0) intensity (see Fig- 
ured]). 

On the other hand, the ratio is not constant in 
the two-phase turbulence model where tempera- 
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Fig. 5. — Ratio of emissivities per unit volume, 
A[cii]/Aco(i-o) as a function of Aco(i-o)- The 
values of A[cii]/Aco(i-o) are shifted by an ar- 
bitrary factor of 3 X 10~^ for clarity of the fig- 
ure. The horizontal line located in the lower half 
consists of dots from the one-phase model, while 
the upper, slightly bending curve consists of those 
from the two-phase model. The thin solid curve 
is the isobaric contour at the saturation pressure 
Psat = 103-4 K cm-3. 



ture strongly varies from place to place. In this 
case, the pressure approaches approximately con- 
stant within the whole calculation domain and the 
points taken from our two-phase model fall on the 
isobaric contour of P = Psat- We overplot the 
isobaric contour with P = Pgat with a solid line 
in Figure El which closely trace the dots from our 
two- phase simulations. Under constant pressure, 
increasing density indicates decreasing tempera- 
ture. In the low density (or high temperature 
T > lOOOK) limit the cooling rate per unit vol- 
ume increases with density as ri^, and the coUi- 
sional cross sections {(jv) are similar for both C+ 
and CO in such a high-temperature regime. Thus 
Aco(i-o) is approximately proportional to Aj^ji] . 
On the other hand, as the temperature decreases 
below T < lOOK, (crw)[cii] decreases with T much 
faster than (o'w)co(J=i) because of the exponen- 
tial decrease in the number of particles that can 
excite the [C II] line (92K). Simultaneous increase 
in density with a decrease in temperature brings 
the number of CO molecules in the level J — 1 
close to the LTE limit. Because of the slower de- 
crease in (o'w)co(J=i) with temperature compared 
with (av) [cii] , and the saturation of the J = 1 level 



population owing to LTE, there is a sharp decline 
of the isobaric contour toward Aco(i-o) ~ 10~^^ 
erg sec~^ cm^-^. This explains the behavior of the 
isobaric line in Figure El which turns downward 
with increasing density. 

3.3.2. CO rotational line ratios 

In this subsection, we consider ratios between 
intensities of C0(1 — 0) and higher J lines. Since 
the excitation energies of rotational levels from 
the ground state Ej/ks = 2.75KJ{J + 1) in- 
crease with J, higher J lines trace higher temper- 
ature gas. In Figure [5] we show the line inten- 
sity ratios Pj,(j_i)/io = ^co(j-(j-i))/-fco(i-o) 
for both two-phase and one-phase models. The 
dependence of line ratios on rotational quantum 
number J is surprisingly different for these two 
turbulent flow models. In the one-phase model. 




10"" 10'" 10"" 10*" 10"^" 
Ico(i-i»)IeTg sec"' cm Str"^] 

Fig. 6. — The ratios of intensities of a CO high 
(J = 2,3,4) excitation line with the C0(1 - 0) 
line. There are two sequences of plots correspond- 
ing to the two-phase and the one-phase model, 
respectively. The green dots are P43/10, the red 
are P32/101 and the black are ^21/10- The left se- 
quence is the two-phase flow data, and the right 
one is the one-phase flow data. In the two-phase 



model, the line ratio R 



J,(J-1)/10 



increases with J, 



while, in the one-phase model, its dependence on 
J is opposite for J < 4. 

the line ratios Pj,(j-i)/io decrease with increas- 
ing J, while this trend reverses in the two-phase 
model up to J < 5: in the one-phase model, al- 
though the CO (2 — 1) line is almost as strong as 
C0(1 — 0) line (P21/10 ^ 1); higher J lines become 
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progressively weaker. On the other hand, in the 
two-phase model, the line ratio -R21/10 is as large 
as 2 and i?32/io and i?43/io become even larger. 

This can be explained by the following argu- 
ment: Since density is lower or at most marginally 
higher than the critical density (rid- ~ SOOOcm^'^) 
in our snapshots, here we consider the low-density 
limit. In other words, we discuss the cases where 
the radiative decay from each level dominates all 
the collisional transitions both into and from this 
level, except the collisional excitation from the 
ground state. Under this assumption, the equa- 
tions of the detailed balance for the optically thin 
lines can be approximately described as 



J>i 



for the ground level (J = 0) and 



(9) 



(10) 



for levels J > 1. The collisional excitation rate 
Cij is written as in equation ([5]). From these equa- 
tions, we obtain 



Aj,j-inj = ^ C( 



0J'"-0- 



(11) 



J'>J 



Equation pip implies that all the CO molecules 
that are coUisionally excited to the state J cascade 
down to lower levels by spontaneous emission with 
A J = 1. Then the downward transition rate from 
the state J by spontaneous emission is equal to 
the sum of the collisional excitation rate from the 
ground state to all the states with J' > J. Us- 
ing equations ([3]) and (fTT|) . the ratio between the 



adjacent lines ( J -f 1) 



i?, 



J and J ^ (J — 1) is 
Ej+i,J Sj'>j+i 7oj 



Ej+i^j Gj+i 



1 + G,7+i ■ 



(12) 



(13) 



E 



where Gj+i is defined as 

7J'0 9j' 

Gj+i = > exp 

7J0 9J 

In the above, we have used the relation 



lOJ 



9J 

7 JO — exp 
50 



knT 



(14) 



(15) 




0.1 1.0 10.0 100.0 1000.0 
nH[cm"^] 

Fig. 7. — Thermal equilibrium curve (F = A, thick 
solid curve), and isothermal lines of T = 5000K 
(dot-dashed), T = BOOK (dashed), and T = 50K 
(dotted : the lowest temperature of the snapshot 
of two-phase model) in the pressure-density dia- 
gram. The saturation pressure P = Pgat is plot- 
ted with the thin solid line. The temperature 
range where bi-stability is possible is limited to 
300K< T <5000K. 



where gj is the statistical weight of the level J . 

In the following, we examine the line intensity 
ratios using the set of equations (fT2|) -(fT4 |) in the 
low and high temperature limits. In the low tem- 
perature limit where T <^ iJj+i^j/fcs, Gj+i can 
be approximated with its first term. 



R 



Ej+i,.J iJ+ifl 9,7+1 
E.i.j-1 "fjfi 9j 



■ exp 



ksT 
(16) 

Equation p6p shows that when the temperature 
is lower than the excitation energy Ej^i_j /ks, 
the line ratio is exponentially re- 

duced. In our one-phase model with T=10K, the 
excitation energy Ej_j-i/kB is higher than the ki- 
netic temperature for high J {J > 2) lines and the 

line ratio Rj,j-i/Rio = Ylj>^i Rj'+ij' 
is significantly smaller than unity. The only ex- 
ception is i?2i/io ~ 1 since Ei,o/kB — UK is 
comparable to the kinetic temperature. On the 
other hand, in the high-temperature limit where 
T ^ Ej^i_j /kB, Gj+i 3> 1 since numerous terms 
in the summation of Gj+i (eq. [14]) are order of 
unity. Then, 



R 



Ej+ij 
Ejj-i ■ 



(17) 
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In oiir two-phase model, high-temperature regions 
where high J Hnes (T > -©21/^3) can be colHsion- 
ally excited are present. In the snapshot shown 
in Figure [U indeed, the temperature ranges from 
50K to 3700K which is higher than CO excita- 
tion energies up to J < 4 {Ej — 55 K) in most 
of the volume (see Fig. [7|). The excitation ener- 
gies up to J < 4 are smaller than 50K. Therefore 
equation (fT7|) is applicable to the whole region and 
R.J. . J ~i/w becomes order of unity. 

3.3.3. [C I] line ratio 

In the transition region between the warm 
atomic region of low density and the cold molec- 
ular region of high density, carbon is partly in a 
neutral atomic form C°. In this section, we calcu- 
late [C I] fine-structure lines. In the temperature 
and density range under consideration, impor- 
tant lines are the forbidden transitions among the 
triplet ^P2, '^^'i and ^Pq in the ground electronic 
state. In Figure [H we show the ratio of ^Pi Pq 
(24K) and ^P2 Pi (39K) hnes. Similar to CO 
rotational line ratios, the [C I] line ratio becomes 
close to unity in the two-phase model, while it is 
significantly smaller than unity in the one-phase 
model. This behavior can be understood in a sim- 
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Fig. 8. — The ratios of C I fine-structure- line in- 
tensities /(3P2 -3 Pi)/J(3pi -3 Pg) as a function 
of /(3Pi —3 Pq). The two sequences correspond 
to the one- and two-phase models. The ratios 
of high- and low-transition-energy lines of C° re- 
semble those of CO molecules: in the two-phase 
model, the ratio is close to unity, while in the one- 
phase model the ratio is smaller than 0.1. See 
discussion in §3.3.21 



ilar way to the argument for CO rotational lines 
( §3.3.2|) . Again, we consider the low-density limit. 
In the following we call levels ^Pq, ^Pi, and 3p2 
levels 0, 1, and 2, respectively. Since the radiative 
transition from level 2 to is negligible compared 
with that from level 2 to 1, the equations for sta- 
tistical balance read 

^2i«2 + Coiuo = Aioui (18) 

for level 1, and 

Co2"-0 = ^21"2 (19) 

for level 2. From these equations and relations ^ 
and (jl5p . the line ratio is given by 



^21/10 — 



-£'2,1^21^^-2 

E2.1 702 
-£^1,0 701 + 702 
G 



= 1.625- 



1 + G 



(20) 
(21) 
(22) 



where 
G 



72032 I E2S 

■ exp 



7io3i 



kBT 
0.12 



0,6 exp(-fi).,23) 



ViooKy 



where we consider only collision with neutral hy- 
drogen as for the excitation partner for numeri- 
cal evaluation. The above expression clearly ex- 
presses that the line ratio becomes close to unity 
for T > 39K (as in the two-phase model), while 
for lower temperature (as in the one-phase model 
of r = lOK) it decreases exponentially. 

3.3.4. Summary of the Line-ratio Analyses 

We have demonstrated in this subsection that 
the intensity ratio Rh/i of lines of high-transition 
energy {AEh/kB « lOOK) and low-transition en- 
ergy {AEi/ks ~ lOK) is as high as unity in the 
two-phase model, while Rh/i is much smaller than 
unity in the lOK isothermal model. This feature 
is quite robust and is independent of the species of 
emitting particles (atomic carbon or CO molecule) 
and the nature of transitions (fine-structure tran- 
sitions or rotational transitions). Discussion in 
^3.3.21 explains the trend in line ratios in terms of 
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the level population in low-density gas. The level 
populations are determined by coUisional excita- 
tion and radiative de-excitation by spontaneous 
emission. In this case the number of particles 
that have enough energy for level excitation is pro- 
portional to exp(— AE'/fcsT), and then the level 
population reflects the ratio of transition energy 
and kinetic energy, AE/ksT. In other words, 
our results show that simultaneous observation of 
both high- and low-energy lines from the diffuse 
atomic/molecular transition region will reveal the 
thermal structures of the gas. If the high line ratio 
is observed in such regions in future, we can con- 
firm the existence of two-phase turbulent flows. In 
the following section, we discuss the implications 
of our results for the CO line ratio in more detail. 

4. Discussion 

4.1. Possible CO abundance variation in 
the virarm phase 

We have assumed uniform abundance distri- 
bution both temporally and spatially for all the 
species. This assumption appears too artificial 
if we consider complex chemical reactions in the 
temperature range T < lO'^K. Actual chemical 
abundances in diffuse clouds will evolve both in 
space and time, reflecting the environments of the 
clouds and their evolutionary history, that is, the 
development of the multi-phase medium. In par- 
ticular, whether CO is present in such a diffuse 
environment is quite uncertain. Since line inten- 
sity is proportional to the abundance of emitting 
species in an optically thin medium, the lack of 
CO in the warm gas would reduce the CO line 
intensity. In Figure [2^, one can see that warm 
(T ~ lOOOK) gas occupies a large fraction of vol- 
ume in the two-phase medium. In evaluating the 
CO line intensity by integrating the emissivity 
along the line of sight (or the x axis in our anal- 
ysis), we need to take into account the possible 
deficiency of CO molecules in the warm phase. 

To assess the extent of this effect on our results, 
we perform the following simple experiments: Al- 
though formation and destruction rates of CO 
molecules are dependent on various thermal fac- 
tors, e.g., density and temperature, the degree of 
depletion onto dust grains and so on in a complex 
fashion, it is plausible to assume that the abun- 
dance of CO, y(CO) increases with decreasing 



temperature and increasing density. In our two- 
phase model, the whole region is approximately 
in pressure balance. Thus higher density means 
lower temperature, and CO abundance can be re- 
garded as a function only of temperature, as long 
as chemical reactions of CO are in equilibrium. As 
the simplest example, we assume that CO is com- 
pletely destroyed above a threshold temperature 
Tcrit. Namely, the CO abundance y(CO) is taken 
as a step function of temperature T; 



y(CO) 



3 X 10 




-7 



T < Tcrit, 
T > Tcrit. 



(24) 



We take Tcrit as a free parameter and examine the 
effect of different values Tcrit on the line ratios. 
Figure [7] shows that the medium is thermally un- 
stable in the range of 300K < T < 5000K. There- 
fore the spread of temperature grows as large as 
this range in our two-phase flow calculations (50K 
< T < 3, TOOK). We study four cases of Tcrit = 
lOOK, 400K, lOOOK, and 3000K, which correspond 
to uh = 40, 10, 4, and 1.3 cm""^, respectively, un- 
der the saturation pressure Tsat — lO'^'^ K cm"'^. 
We repeat the same analyses as in §3 for each of 

Tcrit ■ 

In Figure [HI we show the CO (J — (J — 
1))/C0(1 — 0) line ratios for different values of 
Tcrit- With a decrease in Tcrit, the volume occu- 
pied by the warm phase (where 7;(C0) is set to be 
0) drastically increases. Radiation is emitted only 
from surviving cold clumps. In spite of a drastic 
increase in volume of the warm phase, the line 
ratios do not change significantly (Figure [9]). This 
result is supported by our discussion in H3. 3. 21 that 
the line ratio of CO rotational transitions depends 
only on the ratio of the excitation energy of the 
line AE, and the kinetic temperature of the gas 
/cbT, as long as density is low (n < Jicrit)- As 
depicted by the isothermal lines in Figure [71 after 
the system settles onto a constant pressure around 
the saturation pressure Tsat, the lowest temper- 
ature of the two-phase model is as high as 50K 
in our calculation. This lowest temperature as- 
sures Ej/ksT < 1 for J < 4 (Tj = 55K), so that 
even in the coldest region in our two-phase model 
the line ratios ^j,(j-i)/io exceed (or are close 
to) unity. Since CO molecules are preferentially 
formed in the dense strongly shielded clumps (see 
also discussion in this section below), the temper- 
ature inside the CNM of the two-phase ISM can 
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be lower than our results owing to the resultant 
CO co oling. In fact, according to IWolfire et al 



(|1995[ ). the thermal equilibrium temperature in 
CNM is as low as ^ 20K. The condition for high 
line ratio Ej/ksT < 1 becomes more difficult to 
be satisfied in colder CNM: the line ratio is ac- 
cordingly lowered. Thus, the value of i?j,(,/-i)/io 
of higher J is reduced, while -R21/10 stays the same 
as our results except in strongly shielded molecu- 
lar clouds of lOK owing to the very low excitation 
energy only 16. 5K. In this case we should use the 
line ratio i?j,(j-i)/io of low J lines to study the 
thermal structure of the CNM. 




10" 10" 

Ico(i-o)["S sec ' cm"^ 
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Fig. 9. — CO rotational line ratios with different 
values of threshold temperature Tdit , above which 
CO abundance is set to zero. Dots of each color 
correspond to the same transition J as those dis- 
played in Figure [H] The threshold temperatures 
are (a) Tcrit = lOOK, (b) 400K, (c) lOOOK, and (d) 
3000K. If the pressure is at the saturation value 
-Psat = lO"^'^ K cm~'^, the equivalent threshold den- 
sities are (a) 40 cm"'^, (b) 10 cm^^^, (c)4 cni^-^, and 
(d) 1.4cm~'^, respectively. One can observe that 
as Tcrit increases, the value of /(C0(1 — 0)) in- 
creases. This is because emission from the warm 
gas is added to the summation along the line of 
sight. Also the values of the ratio do not change 
significantly, because, even in the coldest regions, 
the energy difference of the levels AE/kBT < 1 in 
the two-phase model (see discussion in i)3.3.2p . 



The increase in the threshold temperature Tcrit 
is accompanied by an increase in the average tem- 
perature T in the emission region. Since the 
line ratio ^j,(j-i)/io grows larger as the aver- 
age temperature T increases ( ij3.3.2p . the shift of 
the threshold temperature Tcrit does not alter the 
trend toward the averaged i?j_(j_i)/io > 1 as long 
as AE/ksT < 1. Note that our simple experi- 
ments include two extreme cases of high- and low- 
y{CO) values in the warm phase. Thus we con- 
clude that, as long as we are concerned with the 
line-ratio analysis presented in this article, the re- 
sult is not significantly affected by complicated 
chemical evolution. In other words, the line ratio, 
not the intensity of each line, is a good observa- 
tional tool for understanding the thermal struc- 
ture of turbulent interstellar medium. 

Since the chain reaction for formation of CO 
molecules starts with a reaction involving H2, the 
presence of H2 is a necessary condition for CO 
formation. H2 molecules are formed via associa- 
tion of two hydrogen atoms on the surface of dust 
grains. The formation timescale of H2 molecules 
is des cribed as t^.^ = 1.67 x lO^n^^ years tSpitzeil 
19781 ). The formation timescale becomes on the 



order of 10 years in the dense cold clumps, where 
nji ^ 10'^ cm~^. Then if the dense clumps sur- 
vive as long as 10® years, CO molecules can form 
inside. There are some arguments supporting the 
existence of CO molecules in tenuous warm ISM. 
Recently, llnoue. Inutsuka fc Kovamal ( 20061 ) em- 
phasized the importance of gas flows (i.e. evapora- 
tion) from cold phase to warm phase. If molecules 
are abundant inside the cold clouds, they will be 
dispersed into the warm phase. The efficient evap- 
oration of CO from cold clumps enables the supply 
of CO molecules to the warm gas on a relatively 
short timescale, although the chemical timescale 
is much longer than other relevant timescales in 
the warm phase. Of course, if the UV radia- 
tion is as intense as in average diffuse ISM, such 
CO m olecules would be photodissoc iated immedi- 
ately. iHennebelle fc Inutsuka ( 20061 ) proposed an 
alternative model for cold molecular clouds that 
possess warm gas deep inside. In this picture, 
CO molecules can avoid photodissociation within 
such warm neutral medium (WNM) because ex- 
ternal UV radiation is blocked by the surrounding 
cold molecular gas. In this scenario, the condi- 
tion AE/ksT < 1 is satisfied for CO molecules in 
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WNM. 

Here we discuss two possible evolutionary sce- 
narios for CO abundance in diffuse ISM. First, we 
consider the case of relatively strong interstellar 
UV radiation, stronger than the adopted value 
(Go = 1 in our article). If there were a large 
number of FUV photons, CO molecules would be 
easily photodissociated. The shielding of UV pho- 
tons by atomic and molecular hydrogen and dust 
grains would result in selective destruction of CO 
molecules in the tenuous warm medium, but CO 
molecules in dense, cold clumps, surrounded by 
a warm medium, would be protected against UV 
penetration. Another possibility is that a strong 
UV field might therm ally stabilize the m edium at 
a single warm phase (|Wolfire et al. I I1995I) . In this 
case, however, CO molecules would not be able 
to survive photodissociation and CO line emission 
would be very weak, even though there are some 
turbulent flows. The survey observations of molec- 
ular clouds in the Galactic plane suggest that the 
number of regions of high-line r atio i?9i /in > 1. 



is small (I S akamoto et al.l 119971 : 1 Oka et al 



f 1 1991^ . 



These results might be interpreted as a signature 
of the lack of initial CO molecules by photodissoci- 
ation. If this is the case, more refractory molecules 
must be used as tracers of the two-phase medium. 
In order to examine what transitions can be bet- 
ter tracers, it is necessary to include greater detail 
of chemical evolution in our hydrodynamical sim- 
ulations, as well as of the evolution of the main 
cooling and heating sources. Second, we con- 
sider the case of weak UV irradiation as we have 
considered in this article. In our two-phase mod- 
els, the temperatures do not exceed T « 5000K, 
and CO molecules once formed are not subject to 
the coUisional dissociation. In this case, our ex- 
periments in this subsection are straightforwardly 
applicable and the ratio of rotational transitions of 
CO molecules is a good probe of thermal structure 
in the turbulent ISM. 



^Equation (|2]l describes the cooling rate of [C II] 158 /im 
line, the major coolant in tenuous ISM with which we are 
concerned in this article. 



4.2. Comparison with Observational Re- 
sults 

4-2.1. Related Observations oj Diffuse Clouds 

The effective cooling function we adopted 
(equation [21) was calculated under the condi- 
tion of small visual extinction. Ay < 0.1. This 
environment corresponds to the periphery of 
molecular clouds where the composition of car- 
bon changes from atomic to molecular form. 
Sakamoto fc Suna d3 (120031 ) Found peculiar struc- 
tures in the position-velocity (P-V) diagram 
around Ay ^ 1 in their observations of C0(1 — 0) 
in Taurus HCL2, which spans over a wide range 
of Ay. Tiny clumps with sharp velocity gradients 
of radius R « 0.05pc suddenly appeared when 
Ay fell below unity in all of three observed strips. 
They speculated that these structures were cold 
dense clumps formed by thermal instability since 
this instability under mild UV heating occurred 
approximately around Ay ~ 1. In Figure [TOl we 
display the synthetic P-V diagram of C0(1 — 0) 
line from the snapshots of the simulations at 3 
Myr. The two-phase model shows that velocity 
of clumpy structures spreads over a wide range 
(Awr ~ 5 km sec~^). Indeed, our synthetic P-V 
diagram o f the two-phase model looks quite similar 
to that of ISakamoto fc Sunadal |2003h at Ay <1. 
Since our simulations do not consistently include 
the spatial variation of UV shielding in modeling 
periphery of molecular clouds, direct comparison 
with these observational results should be done 
with caution. This similarity is, however, quite 
encouraging for confirming our two-phase model 
driven by UV photoelectric heating. 

IRAM observations of starless molecular clouds 
also suggest the pos sibility that tiny gas c lumps 
are filling the beam (^Fal garone et al. 1998 ). Fur- 
thermore, in more diffuse regions in the Galac- 
tic plane, ISO observations of II2 rotational lines 
suggest the existence of coUisionally heated warm 
gas (T < lO '^K) within the cold diffuse ISM 



(jFalgarone et al.ll2005l ). Although both of the ob- 
served regions are denser than our calculations, 
their structures might be roughly consistent with 
the two-phase mediu m models generated by UV 
photoelectric heating (Ko vama fc Inutsuka 2002t 
lAudit fc Henneb ell c 2005i) or dissipative heatiri g 



of Alfven waves (jllennebelle fc Inutsukal 2006h . 
We can clarify the nature of the multi-phase ISM 
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by making comparisons with line-emission calcu- 
lations like ours and such detailed observations of 
diffuse ISM as described above. 

4-2.2. Comments on observations of the CO line 
ratio i?2i/io 

Despite our conclusion that the two-phase 
medium has a high line-ratio Rjj-i/io = 

^CO(,7-(,7-i))/-^co(i-o) ^ 1-0, observations have 
found that regions of high line-ratio R21/1 



> 



(Sakamoto et al. 1997: Oka et al 



1.0 are of l imite d fraction 
19981 ) and usu- 



ally only appear close to HIT regions. Our results 
of high -R21/10 values are in part due to the as- 
sumption of small optical thickness (r ^ 1). If 
we included the effect of finite optical thickness, 
the line ratio i?2i/io averaged over the simulation 
box would become smaller. However, as for the 
peripheral regions, which are optically thin, our 
treatment is still valid and our conclusions should 
not be altered. Discrepancy between existing ob- 
servations and our prediction is probably due to 
the fact that observations tend to be strongly bi- 
ased to clouds with high CO intensity, thus having 
higher average density compared with the transi- 
tion regions we discussed. To observationally ex- 
amine the atomic/molecular transition regions, we 
must take care not to include both high and low 
CO intensity regions in a single observing beam 
in resolving the spatial gradient of Ay. For this 
purpose, observations not only with high angular 
resolution but also with high sensitivity are needed 
for achieving high S/N observation of faint tran- 
sition regions, avoiding contamination of high in- 
tensity regions. We hope that future instruments 
will provide a good information for the two-phase 
scenario of ISM turbulence. 

Here we discusss some of the characteristics of 
the cold clumps analyzed in this paper and future 
observation. The typical size of the cold clumps is 
determined by the maximum growth wavelength 
in thermally unstable medium Amgr ~ vAfA^xjoI 
(jrield 19651 where Af = (KT/nii^Ay/^ is the 
Field length, and Acooi = c, ■ (3/2)fcBT/(nHA) = 
Cs ■'cool is the cooling length, respectively. Numer- 
ically these values are 
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where nuA is cooling rate per particle, and 
the normalizations are taken from typical val- 
ues in our calculation. If tiny cold clumps re- 
side at the distance of nearby molecular clouds 
D = 140pc, the angular size of a clump of ra- 
dius Amgr = 5.3 x lO^AU is - 3.7". From Fig- 
ures[6l[8]and equation ([7]) intensities of C0(1 — 0) 
and [C I] lines of the two-phase medium are 
/(min) > 10~^^(Aa;/2.4pc) erg sec~^ cm~^ str~^ 
in our calculations (where Aa; is the total length 
scale of the two-phase medium along the line 
of sight). This intensity value corresponds to 
the brightness temperature Tf, = 2.46 x 10~'*K 
(Aa;/2.4pc)K, and is too weak to be detected with 
ALMA telescope unless the line-of-sight length of 
the cloud is sufficiently large (ALMA Sensitivity 
Calculatorjl, even though its high angular reso- 
lution is sufficient for resolving each tiny clump. 
This might suggest that we should target near-by 
cloud s such as some of the hig h latitude clouds 



fe.g.. iMagnani. Blitz fc Mundv 1 11985). In addi- 



tion to the weak average intensity, the volume fill- 
ing factor of the tiny CNM clumps of which inten- 
sity dominates, is much less than unity (Fig. 
Taking into account of these difficulties, we should 
be careful in the interpretation of the observations 
for the tiny CNM clumps in the diffuse two-phase 
ISM even with ALMA. 

On the other hand, Figure [TO] shows that 
the tiny clump-like structures are assembled into 
larger structures of < 0.1 pc when integrated along 
a line of sight in the P-V diagram. The flux den- 
sity of such assembly is 
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Fig. 10. — Synthetic position-velocity diagram of 
the C0(1 — 0) Une at ^ 3Myr snapshots of the 
two-phase model. Tiny structures spread with a 
velocity dispersion of Avr ~ 5 km/s. 



which might be observable with a future telescope 
equipped with highly sensitive receivers. For bet- 
ter evaluation of the intensity, we need to study 
the global structure of the atomic/molecular tran- 
sition region which supposedly consists of the tiny 
clump- like structures appeared in our simulations. 
Such a global modelling is beyond the scope of this 
paper, but should be included in future simula- 
tions. 

When future highly sensitive multi- J line obser- 
vations of "peculiar" structures with typical size 
O.lpc detected in the strip-s canned C0(1 — 0) data 
( Sakamoto fc Suna dal 120031: see §4.2.1) are avail- 
able, they will greatly improve our understand- 
ings of their origin and identity. In those rims 
of molecular clouds with Ay < 1, even smaller 
structures might be found with future deep obser- 
vations. Small (< 0.1 pc) structures with high- 
line ratio in imaging observations would be the 
best verification of the two-phase ISM turbulence 
models and of the formation scenarios of molecular 
clouds b ased on the multi- phase turbul ence mod- 



els (e .g., iKovama fc Inutsuka 2000: Bergin et al 
20041 ). Detailed observations using forthcoming 



powerful instruments will reveal the thermal struc- 
ture in diffuse interstellar medium, as well as in 
dense molecular clouds. We hope our calculation 
of the carbon line ratio would be a prototype to 
study the nature of ISM using forthcoming high 



quality observational data. 

5. Summary and Conclusion 

We examined observational characteristics of 
the two-phase turbulence model by combining hy- 
drodynamical simulations and radiation calcula- 
tions. We modeled the diffuse interstellar medium 
irradiated by UV radiation of Gq = 1, attenuated 
with an average visual extinction of Ay ~ 0.1. 
In our model settings, the two-phase structures 
are formed by photoelectric heating, and therefore 
our model corresponds to the atomic/molecular 
transition region at the periphery of dense molec- 
ular clouds. In order to examine such transi- 
tion regions, we calculated line intensities based 
on nonLTE level populations for three species of 
carbon, C+, C°, and CO. By comparing two line 
intensities of high- and low- temperature tracers, 
we found that the line ratio Rh/i distinctly re- 
flects the thermal structure of the medium. The 
results of our analyses showed that, in the two- 
phase model, higher transition energy lines can 
be stronger than those of lower energy as long as 
AE/ksT < 1. Our argument on the level popula- 
tion in low-density limit (n < ricrit) demonstrated 
that the line ratio of high- and low- temperature 
tracers are determined by the ratio of transition 
energy and gas kinetic energy, AE/ksT. In ad- 
dition to uniform abundance calculations, we also 
analyzed cases where CO is dissociated above a 
threshold temperature, or equivalently, below a 
threshold density. Our analyses showed that abun- 
dance distribution of CO molecules in the warm 
phase does not significantly affect the high line ra- 
tio that characterize the two-phase model, as long 
as CO molecules are not completely photodisso- 
ciated. This is because, even in the coldest re- 
gions in our two-phase model, the temperature 
is as high as 50K (« E^/kB, the energy level of 
J = 4 of the CO rotational transition), and sat- 
isfies AE/ksT < 1. By seeking these features 
in the diffuse atomic/molecular transition regions, 
we will be able to confirm the existence of two- 
phase interstellar turbulence by using the high an- 
gular resolution and high sensitivity of future in- 
struments. 
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